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Extra-galactic sources of photons have been used to constrain space-time quantum fluctuations in the Uni¬ 
verse. In these proposals, the fundamental “fuzziness” of distance caused by space-time quantum fluctuations 
has been directly identified with fluctuations in optical paths. Phase-front corrugations deduced from these 
optical-path fluctuations are then applied to light from extra-galactic point sources, and used to constrain var¬ 
ious models of quantum gravity. However, when a photon propagates in three spatial dimensions, it does not 
follow a specific ray, but rather samples a finite, three-dimensional region around that ray — thereby averaging 
over space-time quantum fluctuations all through that region. We use a simple, random-walk type model to 
demonstrate that, once the appropriate wave optics is applied, the averaging of neighboring space-time fluc¬ 
tuations will cause much less distortion to the phase front. In our model, the extra suppression factor due to 
diffraction is the wave length in units of the Planck length, which is at least 10^^ for astronomical observations. 


PACS numbers: 95.75.Kz, 03.67.Lx, 04.60.-m, 98.62.Gq 

I. INTRODUCTION 

There have been several proposals to use extra-galactic 
point sources to constrain the quantum fluctuations in space- 
time EH. It was argued that, space-time fluctuations cause 
random phase shifts of photons, and that these shifts accumu¬ 
late throughout the very long light propagation path from the 
point source to the earth, causing wavefront distortion from a 
perfect spherical shape upon arrival at the earth. The manner 
in which these fluctuations accumulate depends on the spe¬ 
cific model of quantum gravity phenomenology; and in par¬ 
ticular, it has been claimed that the random walk model could 
be ruled out by existing imaging data from the Hubble Space 
Telescope. In this paper, we point out a serious omission in 
the theory so far employed by all such proposals, and argue 
that the random walk model, once given a closer look, cannot 
be ruled out at all by current or any forseeable observations of 
extra galactic sources of photons. 

In Refs, nail, based on the argument about the quantum 
space-time on the phenomenological level luia, the authors 
assumed photons originating from a point source to undergo a 
random phase shift Acj) due to space-time fluctuations: 

A(/) ~27r(/p/A)“(L/A)'^“, (1) 

where (p is the Planck length, A is the wavelength of the 
light, L is the light propagation length, and a is a parame¬ 
ter that depends on specific models of quantum space-time 
phenomenology. In particular, a = 1/2 corresponds to the so- 
called random-walk model, which can be understood as hav¬ 
ing the speed of light fluctuating dramatically at the scale of 
/p. The relation ([T]) stems from phenomenological description 
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of the quantumness of space-time by assuming that the uncer¬ 
tainty in distance measurement 5L due to space-time fluctua¬ 
tion over a distance L is given by 

5L > , (2) 

which was first derived by Ng and van Dam and then dis¬ 
cussed by others lEHll. Note that Eq. Q is related to Eq. E 
byA0 -27r(5L/A). 

If we consider a photon propagating along one spatial di¬ 
mension. Suppose we divide its propagation distance L into 
pieces of /p; within each piece, the phase-shift fluctuation it 
gains is substantial; 5^ ^ 2nl^/X —while fluctuations in dif¬ 
ferent intervals are independent of each other. In this way, the 
total phase shift of the photon does a “random walk” while 
the light propagates. At the end of propagation, we have a 
photon-phase fluctuation of 

(A0)iD~V^5(/)~-y/pL/A2, (3) 

which is exactly Eq. Q with a = 1/2. 

However, controversies about the magnitude of this wave- 
front distortion effect were raised in several papers. Coule 
Cl first qualitatively argued that this decoherence effect 
would be very small because of the van Cittert-Zernike the¬ 
orem 0. Later, by modeling the discrete space-time in a 
Lorentz-invariance way, Dowker et.al OD calculated a sim¬ 
ple model describing the electromagnetic potential generated 
by oscillating charges at the source rather than studying the in¬ 
dependent dynamics of the light wave field. They found that 
the signal due to space-time foam would be undetectably tiny. 

In our paper, we focus on the independent dynamics of 
light wave fields in space-time with Planck-scale fluctuations 
by quantitatively studying a model. Actually the light field, 
as a wave, does not propagate in the way suggested in m, 
but in the following Huygens-Eresnel way: when a photon 
travels through a space-time region, it does not follow only 
one particular ray, whose length might be subject to the fun¬ 
damental “fuzziness” prescribed by Eq. E, but instead, the 
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wave nature of light, or the quantum mechanical nature of the 
photon, dictates that the photon would simultaneously sam¬ 
ple an ensemble of many different neighboring rays, each of 
which has a potentially different realization of the fundamen¬ 
tal length fluctuation; the actual path-length fluctuation must 
then be given by an averaging among these different length 
fluctuations. This allows to go below IndL/X. Moreover, 
because the linear size of the sampling region can be much 
bigger than Ip (the correlation length of fundamental quantum 
fluctuations), this averaging can dramatically suppress the ac¬ 
tual from 2tz5L/X, or Eq. Q. 

In other (simpler) words: (i) propagation of photons is de¬ 
scribed by the photons’ wavefunction, which describes elec¬ 
tromagnetic (EM) waves; (ii) diffraction of EM waves makes 
them insensitive to fluctuations and disturbances at scales 
much less than the wavelength; and (iii) as an EM wave prop¬ 
agates through a large distance, the coherence level of its 
phase front increases unless further disturbance keeps com¬ 
ing ||9| (This point was also independently discussed quali¬ 
tatively by D. H. Coule ifTOlH . If we divide 3-D space into 
cubes of side length ~ Ip, with disturbance to light propa¬ 
gation independent within each cube, then using Eourier op¬ 
tics, it is easy to estimate that the effect we have described 
above will be suppressed by yJlp/X in each transverse direc¬ 
tion (because only perturbations with spatial frequencies be¬ 
low ^ \/X along each transverse direction get registered by 
the propagating light), which leads to 


quantum fluctuations. In addition, light-speed fluctuations in 
regions separated by more than the Planck length are inde¬ 
pendent of each other following the random walk model. As 
we shall see in the calculation below, the particular shape of 
the correlation function Eq. (|^ does not matter — the effect 
will remain the same as long as: (i) the total variance in n 
is of order unity, and (ii) coherence in e exists only between 
points separated by less than than fp. (3) Our toy model entails 
a quantum gravity induced breaking of Lorentz invariance 
which has been strongly restricted by astrophysical data HI. 
However, our toy model should be sufficient to demonstrate 
the omitted effect in CHS due to the wave nature of light and 
to capture the key character of the previously proposed idea. 


III. THE WAVE EQUATION 

Returning to the wave picture, we first write down the wave 
equation: 

-[l+2e(v,y,z)]^+v24> = 0. (7) 

Since our refractive-index perturbation is static, we can ex¬ 
pand the total wave into two monochromatic pieces, the ideal 
wave fPo{x,y,z)e^‘“°‘ and the scattered wave \j/{x,y,z)e^'“°‘: 

^=[^o{x,y,z) + w(x,y,z)]e-“^\ ( 8 ) 


A03D = A(/)id(/p/1). (4) 

In the following, before discussing the consequences of this 
rather dramatic suppression, we shall justify Eq. 0 at a very 
pedagogical level. 


II. MODEL OE SPACE-TIME ELUCTUATION 

Let us construct a toy model reflecting the effect of quan¬ 
tum space-time fluctuation induced length uncertainty on the 
propagation of light. In this model, the Minkowski space- 
time with length fluctuation at each point is viewed as a 
“medium” with a random (yet static) distribution of refrac¬ 
tive index n{x,y,z) = 1 + e(x,y,z). We assume the following 
translational invariant spatial auto-correlation function for e: 

{e{x,y ,z)e{x ,y ,z )) 

= a^n{x'-x")n{y'-y")n{z:-z"), (5) 

where “(...)” stands for ensemble average, a is of order unity, 
and 


The ideal wave is the unperturbed part of the light field sat¬ 
isfying: —(9^0o/(9f^ -I- V^4>o = 0. At leading order in £, we 
have 

(V^ -p col) W(x,y,z) = -2(ols{x,y,z)^o(x,y,z) , (9) 

which means y/ is a perturbative field sourced by a beat be¬ 
tween the ideal wave and space-time perturbations. We know 
from the observed phenomenology that such a perturbation 
must apply to our situation: namely, fluctuation caused by 
the space-time foam is indeed very weak compared with an 
ideal wave propagating across the Universe, and we have 
\\I//^q \< 1 . 

Eor a point source, we assume 

g'®0'' _ 

<Po{x,y,z) = -—, r=ffx^+y'^+z^. (10) 

47tr 

At the distance L, the scattered wave \j/ must be compared 
with the ideal wave to characterize the modulation caused to 
the idea spherical wave by space-time perturbations. Let us 
define 

a + i(j) = ^ =4-nLYe^‘“'’^, a,^ gR, (11) 

<Po 


n(x) 


1 |x| < Ip 
0 jxj > Ip. 


( 6 ) 


Our purpose is to study the propagation of light as a scalar 
wave traveling in this medium. In this specific toy model, the 
coordinate speed of light propagation fluctuates, in on a small 
region with the size comparable to the Planck scale — this 
simulates light propagation in a space-time with Planck-scale 


so a describes the amplitude modulation, and <j) describes 
phase modulation in radians. We also define the total mod¬ 
ulation, 

^ = ^{a^ + 02 ) = , ( 12 ) 

whose standard deviation is greater than those of both the am¬ 
plitude and the phase modulations. 
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IV. SUMMING OVER PATHS 


We have therefore confirmed Eq. Q. 


To arrive at the answer quickly, we use the Huygens- 
Fresnel-Kirchhoff scalar diffraction theory, which is equiva¬ 
lent to applying the outgoing Green Function, and obtain 112: 


V/(x) = f -2a^e(x')<I>o(x') 

J|x'|<L 


giojolx-x 


4;r|x —x' 


■dx' 


/ 

J\x'\<L 




47i:|x' 


-[-2®o^e(x')] 


gjfllolx-x' 


47i:|x —X' 


-dx'. (13) 


V. SPATIAL-SCALE CUT-OLE 

To separately study fluctuations at different spatial scales, 
we solve the same problem by decomposing the secondary 
wave into modes: 

\l/(r,9,(j)) = Y,l\l/em(i-)Yemi9,(p)]- (19) 

£m 


Note that we have considered only contributions from fluctua¬ 
tions at distances smaller than L to the point source. The inte¬ 
gral ( [T3 ] i can be interpreted as a path integral — over all paths 
that consist of two straight sections (each associated with a 
propagator), and a deflection in the middle due to interaction 
with refractive-index fluctuations (associated with a coupling 
coefficient). Paths with more than one deflection do not have 
to be taken into account in our perturbative treatment at lead¬ 
ing order. 

If we discretize the integration domain, a sphere with vol¬ 
ume ~ L^, into cells with linear size ^ Ip and volume vp ^ Ip, 
we will get a total of Ntot /^p individual cells, each of 
which has a statistically independent fluctuation in e with vari¬ 
ance [cf. Eq. Q]. Then the fluctuation variance of \j/ given 
in Eq. ([T3ll can be estimated in the following way; 




dxdx! 


g;(i)D(|x'| + |x-x'|-|x"|-|x-x"|) 

|x-x"||x-x'||x'||x"| 


(e(x')e(x")), 

(14) 


where the integration is over the region |x| <L,|x'| < L. Since 
the correlation function of £(x) only contributes to the integral 
when |x'-x"| < \lp\, therefore the exponential in the numer¬ 
ator is approximately equal to 1. If then we change the inte¬ 
gration argument using Ax = x' — x" and 2xo = x' -f x" while 
substituting Eq. Q, the variance becomes: 


{W? 


®o 


c/^Ax 


|Ax|<//> 


'|xo|<L 


d^xpj 


Xo|2|x-Xo|2' 


(15) 


The first integral on Ax is just Ip while the second one on xq 
can be approximated to be /t!^. Using the definition of vp 
and Ntotah finally we have: 

2 

~ (16) 

According to Eq. ( [T^ and comparing with Eq. (|^, we have 
A(/)<^~flCUo2y^-(A(/))iD(/pA). (17) 


There is an extra suppression factor of Ip /A, which arises from 
the fact that in Eq. ( [T3] l, the intermediate point x' of the optical 
path has the freedom to depart away from the axis connecting 
the source point and the field point, and sample through Atot 
A//p independent fluctuations, instead of only Np ^ L/lp in 
the one-dimensional treatment. A more precise calculation 
gives only an additional numerical factor of the order of unity: 


A(j) 




(18) 


Here ¥£^{9, <p) are spherical harmonics. They describe angu¬ 
lar variations at scales of Injt, at a radius r, that corresponds 
to transverse length scales of 2nrll, or transverse spatial fre¬ 
quencies of Iji^nr). The modal decomposition of Eq. is 


r dr 



m- 


£(£+ 1 ) 




Inr 


( 20 ) 


with 


r27l pK 

£im{r)= d(p sin9de[£{r,e,(p)Yl^{e,(p)] , (21) 

Jo Jo 


which satisfies 


{£emir)ep^,{r')) = r'j/p/r^ . (22) 

Here we simply assumed 

(e(x)e(x')) = A/p5(^^(x —x'). (23) 


The spatial spectrum corresponding to this correlation func¬ 
tion is identical to that in Eq. © at low spatial frequencies, 
but continues to exist in orders higher than 1 /Ip. In principle, 
those modes may also add incoherently to our output fluctua¬ 
tions, but as we shall see, their contributions will be negligi¬ 
ble. 

Solving Eq. ( [20l i, assuming regularity at r = 0 and the out¬ 
going wave condition at r = L, we obtain; 

3 (0 

Vem(L) = ['^dr[rMcOor)e"^'-eeM 

2n Jo 

(24) 

where je and are spherical Bessel and first-kind spherical 
Hankel functions 1201 . From Eqs. ( |22] i and ( |24l i, we obtain 

leK^L^iy/MWkt') = , (25) 


with 


2 pL 

^|„ = 4A COoLh^PiL) {(Oolpf jj{(Oor)cOodr. (26) 


Note that is independent of m, which is a consequence 
of the rotation invariance of the refractive-index fluctuations. 
The total fluctuation sir = L will then be 






E 


2£+l 

47t 




(27) 
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FIG. 1: Plots of as functions of £/{coqL), for cases with coqL =100 
(dotted curve), 200 (dashed curve), and 300 (solid curve). 


Physically, describes fluctuations at the angular scale ^ 
Inji, or transverse spatial scales InLjl, or transverse spatial 
frequency of il{2nL). Inserting Eq. (|2^, we have 


= a^{(OQhf 

2£+l 


(OqUi'P {(OqL) 


fOJoL 


jl{R)dR. (28) 


It might seem difficult to evaluate the summation ( [28] l ana¬ 
lytically. But since we are solving exactly the same problem 
as the previous section, it should be clear that [cf. Eq. ®] 





(A(/))id(/pA) , 


(30) 


as we have verified numerically in the cases (OqL = 100, 200, 
and 300. 


VI. DISCUSSIONS 

So far in this paper, we have calculated fluctuations on 
the phase front of an extra-galactic point source, caused by 
Planck-scale fluctuations in refractive index — a toy model 
for space-time foam. If diffraction of light were ignored, or if 
we assumed a space-time with one time dimension plus one 
spatial dimension, our toy model would give comparable re¬ 
sults to previous esimates on the random-walk model claimed 
in EHS- However, the diffraction of light requires us to aver¬ 
age space-time fluctuations over all different possible optical 
paths that extend to all three spatial dimensions. This averag¬ 
ing Alters out all fluctuations with transverse scales finer than 
the wavelength. In our model, this causes an extra suppression 
factor of /p/A, with [cf Eq. ([T])] 

A0 < y(pL/A2(/p/A). (31) 


We expect f ^ (OqL, or 2;rL/£ ~ A to be the turning point, 
because at this point the transverse spatial scale is comparable 
to the wavelength A. 

Mathematically, for £ < (OqL, the spherical Bessel and Han- 
kel functions are wavelike at r ~ L, indicating propagating 
waves', for £ > (OqL, the spherical Bessel and Hankel functions 
are not wavelike at r ^ L, indicating evanescent waves. In the 
limiting regimes of ^ <C (OqL and £ ^ (OqL, can be evaluated 
analytically, using asymptotic expansions of spherical Bessel 
functions; 

_ f 1/2, f<a)oA, 

a^imlpy \ (OoL/l(2£+lf7(], £:^(OoL. 

Note that not only does approach 0 at orders £ ^ (OqL, 
the summation of all these higher modes also gives a negligi¬ 
ble contribution. This qualitatively confirms a cut-off at the 
transverse scale of A: fluctuations at much finer scales do not 
generate a secondary wave. This justifies our original use of 
Eq. ( [2^ , and also ensures that the detailed shape of the cor¬ 
relation function (|^ does not matter. In Eig. [T] we study the 
transition zone of £ ~ (OqL numerically, for moderately large 
values of (OqL = 100, 200, and 300, by plotting as a func¬ 
tion of £/ {(OqL). As (OqL -1-00, asymptotes to a smooth, 
non-zero function for £/{( 0 qL) < 1, and to 0 for £/ {(OqL) > 1. 
This means, in the realistic situation of (OqL <^1,£/{(OqL) — 1 
is a sharp turning point between propagating and evanescent 


With respect to previous literature, this suppression is at least 
by 29 orders of magnitude for astronomical observations, if 
the wavelength of A = 10^® m were to be used. Numerically, 
we have 

A(/) ~ 10^26^^/Qp^ ("5 X if , (32) 

which makes the quantum foam effect on the light propaga¬ 
tion extremely small. Suppose we consider another type of 
experiment, namely the use of y-ray arrival time. It is straight¬ 
forward to convert 

M=^<^-=tp^JlvL/l^ = \0-^^s^jL/Gpc{Ey/GoN). 

(33) 

Eor both methods, the effect of random-walk-type space-time 
foam is too small to be seen. 

Eor certain extra-dimension models d, e.g, the one raised 
by Arkani-Hamed et. al d, in which the fundamental scale 
of nature is the electro-weak scale Iem 10^'^m, Ip is am¬ 
plified by a large factor, therefore apparently increasing the 
detectability of quantum foams. Using Ip —)■ Iem ~ 10^^^m, 
we update the above estimates to 

A(j)ed ~0.1 VA/Gpc[(5 X 10-"m)/A]2, (34) 

and 


waves. 


Ated ~ (10^^ s) ^JL/Gpc{Ey/GoW). 


(35) 
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At first sight, these seem more promising to detect. How¬ 
ever, we must be more careful in connecting A 0 ed and Afgd, 
which are the total variance after summing over all spatial fre¬ 
quencies, to observables in actual experiments. In both cases, 
from our study of wave propagation, we have 

(36) 

I 0 |x —x'l > A, 


and a similar relation for the two-point correlation function of 
At{x). 

For a telescope image, the ideal wave gives an Airy pattern 
on the focal plane, while the scattered wave would create a 
diffuse background on the focal plane, which has a total en¬ 
ergy of A 02 ^, and therefore a flux kg^und *at is 


background 




image 



(37) 


compared with the typical flux of the ideal image. It could 
be almost impossible to detect such a background due to con¬ 
fusion with other types of background. For example, we can 
estimate the magnitude of ( [37] l for a typical quasar source 
PG2112H-059, which is ^ 2.5 Gpc away from us. Accord¬ 
ing to the observational data collected by using the Hub¬ 
ble Space Telescope (Wide Field Camera 3), at wavelength 
A ^ 1.25 X 10^® m, the ratio between photon flux of the sky 
background and photon flux of the image is ~ 0.7 x 10 
However, ^b^ackground/'^image is ~ 0.8 X 10^'^, which is 10*^ 
smaller. 

For detecting arrival time of 7 -ray photons, we have to be 
aware that the true Af also depends on the averaging area of 
our detector — which corresponds to the pixel size, which we 
also denote by . As a consequence, we have 


^Gbs — ^^ed 



(38) 


which is likely to gain an additional factor from Afgd. The 
actual time resolution of a gamma ray detector, e.g. Fermi 
Telescope, is ^ 10 /is HD. From (36), the observed correction 
to the arrival time of gamma ray photons with Ey=\ GeV and 
source distance ^ 1 Gpc is ^ 10^^^ s, which is also too small 
to resolve. 

What has been left out in this paper is possible fluctuation 
in time. If the refractive index fluctuation has a white noise 
spectrum up to the Planck frequency, then by observing a por¬ 
tion of the optical spectrum A®, the effect will be suppressed 
further by a factor of 


\/fpA®, 


(39) 


which is likely to be very large as well. 
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